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Abstract 



o 



In this work, we consider frustrated network of cosmic strings to explain possible deviation 
from ACDM behaviour. We use different observational data to put constraint on the model and 
show that a small but non zero contribution from the string network is allowed which can explain 
the possible small departure from ACDM evolution. By calculating the Bayesian Evidence, we 



o ■ show that our model and the concordance ACDM model are equally favored by the observational 
data. 
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1. Introduction 

>. 

tJ" . One of the greatest mysteries in cosmology today is the nature of dark energy in the Universe. 

^ i According to the concordance cosmology, this makes up around 70% of the total energy density 
"Tj- ; of the universe and causes the universe to accelerate around the present epoch Yet, there 
■ is no strong theoretical as well as observational indication to pin point the nature of this dark 
^ ' component. In fact we still do not know whether the mysterious late time acceleration is due to 
y—i . the presence of dark energy or due to any modification of gravitational laws on large scales. 
j>I ! Nevertheless, majority of the observational data j2|, [s], Hf, S] support the concordance ACDM 
model as a possible explanation for the late time acceleration. In this model, the presence of 
the cosmological constant (A) poses two serious problems. First one is related to the fine tuning 



5^ ■ of the value of A to a ridiculously small number (~ 10~^^°Mpj in natural units) necessary for 
the present day cosmic acceleration. The other one is related to its dominance over the matter 
component of the universe precisely at the present epoch. Any small departure from these two 
conditions will result a cosmological evolution that is different from the observable Universe. 
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These conceptual problems dilute the superiority, the ACDM model enjoys over other dark 
energy models as far as observational results are concerned. 

Interestingly, current observational data also allow deviation from the ACDM behaviour. 
Although the deviation is not large, but it is still detectable with the current and future precisions 
of the observational set ups. Hence efforts are on to construct models that can explain this 
deviation. Almost all the approaches for this purpose, assume that A is exactly zero in our 
universe, and an evolving dark energy is solely responsible for the late time acceleration of the 
universe. Possibilities for this include quintessence jsj, k-essence an arbitrary barotropic 
fluid like Chaplygin gas or its various generalization and many more jsj. All these models are 
generally phenomenological without any strong theoretical motivations and they also burden us 
with the challenge of explaining several other model parameters. 

The other option is to assume the existence of a small but non-zero A as in ACDM model 
but to consider some extra component to explain the observed deviation from ACDM behaviour 
joj. This is motivated by the fact that there is no mechanism available that can make A vanish 



without any conflict with the meaningful physics [10[. Moreover the string-landscape model may 



explain the existence of a small but non zero A in our present universe . But it is important 



to have well motivated candidate that can explain the deviation from ACDM; otherwise this 
approch will have the same shortcomings that are present in standard dark energy models. 

In this study, we propose for such a model. There are models where dark energy is assumed 
to be a perfect fluid having equation of state w = ^ < 0. But one major problem for such models 
is due to its sound velocity = w < which causes instabilities on the^small scales. It was 
then proposed to include rigidity in the fluid to make it an elastic solid [12]. In this case, if the 
rigidity is sufficiently large, cj. > even if w < 0. 

A solid perfect fluid can arise in different theories. One such posibility is the frustrated 



network of cosmic strings [13|, Il5|, [IGj. These are topological defects formed during the phase 
transitions in the early universe and can survive till present day. If these strings are produced 
during the phase transition at the electroweak scale, they are very light. But there can be many 
such strings and their combined energy density can be finite. If the cosmic strings intercommutes, 
their evolution obeys a scaling solution. In that case, the energy density scales as in the 
matter dominated era or in the radiation dominated era. In such scenario, strings can never 
accelerate the universe. On contrary, if the strings do not intercommute or pass through each 
other, then the network can be frozen in comoving coordinates. These are called frustrated 
network of strings. Such network of strings will have the equation of state w = —1/3 and hence 
their energy density will go as ps ~ a~^. A phase transition around few TeV and a string 
separation today of the order of few A.U. will result ilg today in the range of and l[l3j . So 
it is possible to get a universe containing fraustrated network of cosmic strings that can explain 
the possible small deviation from the ACDM. For this we do not need Qs today to dominate 
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the energy budget of the universe; a small but non zero value can be sufficient for this purpose. 
This is certainly achievable from a well motivated particle physics set up. Also note that these 
network of strings are not theoretical fantasies; they can be seen in the biaxial nematic liquid 
crystals [17i]. Obviously the question is whether in any reasonable particle physics model, we 
get a fruastrated network of strings. It was shown by several numerical simulations that abelian 
string networks intercommute and give rise to scaling solution without frustration 18|| . But later 



on it was shown that in non-abelian string network, intercommutation can be prevented due to 
topological constraints and they can form stable frustrated network [13']. 

There is also another simple way that can result a string-like term {ps ~ a~^) in the Hubble 
equation. A fluid with a constant negative pressure ( p = —po, with po being a positive constant) 
behaves like a ACDM model. This can be easily seen by putting the constant p in the energy 
conservation equation p + 3H{p + p) = and then integrate it out to obtain the expression for p. 
Let us model small deviation from the ACDM behaviour around present day [z = 0) by Taylor 
expanding around this constant p behaviour: 

P = -pO + PlZ + P2Z^ , (1) 

where we keep terms upto second order, pi and p2 are related to ^1^=0 and ^^\z=o respectively. 
With this one can integrate the energy conservation equation to get the expression for energy 
density and pressure as a function of scale factor: 

/ N / , ^ c 3 pi - 2p2 3p2 
p{a) = {po+Pi-p2) + ^-- Y 

p[a) = -{po+pi-p2) + + (2) 

Here c is an arbitrary integration constant. For the choice pi = 2p2 and p2 < 0, these gives a 
mixture of fluids containing A, matter and a fraustrated network of cosmic strings. 

Hereafter we consider such a mixture of fluids and discuss its various cosmological implica- 
tions. We start with the Einstein equations: 



-3H^ ^^],oa-'('+"'^) - ^a] (3) 

Here subscript "0" means values at present {z = 0, qq = 1), Wg is the equation of state for the 
network of strings. Flatness condition demands that Qmo + + ^so = 1; hence two of these 
parameters are independent. Henceforth we shall consider Q^q and Qgo- As we have discussed 
in the previous paragraphs, for a frustrated network of cosmic strings Wg = —1/3. Usually the 
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equation of state of a cosmic string network is related to its root-mean squared (rms) velocity 
as 14 



ws = \{2vl~l). (4) 

In the limit Vs — )■ 0, one gets the equation of state for string network as = — 1/3. It 
has been shown through different numerical simulations that in the string dominated universe 



Vg < 0.17 [19[ but can never be exactly zero for nonintercommutating strings. In this case, it 
can be shown that the string network can never be conformally streched to the Hubble scale 
preventing them to fraustrate. But in our set up, the energy density of the string network is sub- 
dominant and never governs the expansion of the universe. Hence even if ^ 0, the correlation 
length can grow much slower than the expansion rate of the universe and the network can be 
conformally streched to Hubble scale and can fraustrate. This possibility was also discussed by 
Sousa and Avelino |20|. In our subsequent study, we assume the parameter v"^ to be non zero. 



and see how far one can constrain this parameter together with the other model parameters like 
Wg, ^mo and Qsq using the presently available observational data. 

The equation of state for the combined dark energy fluid ( A+Cosmic String) is now given 

by 

_ Qmo + ^^so - 1 + l{2v^s - l)^^.oa-^(^+"'^ 

Later, we reconstruct this dark energy equation of state, Wde, to see how much deviation from 
w = — 1 is allowed by the current cosmological observations. 

As we discussed in previous paragraphs, if we consider this network of strings as an elastic 
fluid by introducing rigidity, one can get > for this elastic fluid which will ensure stability 



in the small scale perturbations. It has been shown by Battye et al. |2l| that under reasonable 
but not totally general circumstances, one can obtain = 1/5 for a elastic network of cosmic 
strings. In our subsequent calculations, we shall assume this value for for the cosmic strings. 

With this dark energy set up, we use the present observational data to constrain our model. 
We begin by considering the Type la Supernovae observation which probes directly the cosmolog- 
ical expansion. These observations measure the apparent luminosity of the Supernova explosion 
as observed by an observer on earth. In cosmological terms, this is given by luminosity distance 
d^i^z) defined as 
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Using this, we evaluate the distance modulus /i (which is an observable quantity) as 

;, = m-M = 51og-^ + 25, (7) 
Mpc 

where m and M are the apparent and absolute magnitudes of the Supernovae respectively. We 
consider the latest Union2.1 data compilation consisting of 580 data points for the observable 

Next, we use the observational constraints on Hubble parameter as available in the literature. 



Recently Moresco et al. [22| have compiled the Hubble parameter measurements in the redshift 
range < z < 1.75 using the differential evolution of the cosmic chronometers. They have built 
a sample of 18 observational data point for H{z) spanning almost 10 Gyr of cosmic evolution. 
These values are given in Table 1. To complete the data set, we also use the latest and the most 
precise measurement of the Hubble constant Hq j^. 



H{z) (Jh{z) Ref. 



0.090 69 12 [23 
0.170 83 8 [23 



0.179 75 4 [22 
0.199 75 5 [22 



0.270 77 14 [23 
0.352 83 14 [22 



0.400 95 17 [23] 

0.480 97 62 [2^ 

0.593 104 13 [22] 

0.680 92 8 [221 



0.781 105 12 [22 
0.875 125 17 [22 



0.880 90 40 ]24\ 
1.037 154 20 [22.1 



1.300 168 17 [23] 

1.430 177 18 [23] 

1.530 140 14 [23] 

1.750 202 40 [231 



Table 1: H{z) measurements (in units [kms ^Mpc ^]) and their errors. 



Next, we use the combined BAO/CMB constraints as derived recently by Giostri et al. 28 
For this, one starts defining the comoving sound horizon at the decoupling as: 
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1/(1+^.) 



da 



m{a)^TTWh^l^l^'' 



(8) 



where VL^q and ^2^0 are the photon and baryon density parameter respectively at present [z = 0) 
and we have assumed c = 1. Here is the redshift at decoup ling and is accurately given by 
the formula as obtained in 29|. In accordance with WMAP7 j26[], we put = 1091 exactly. 
Another important quantity is the redshift of the drag epoch [z^i ^ 1020), when the photon 
pressure is no longer able to avoid gravitational instability of the baryons. 
Next we define the "acoustic scale": 



TT 



dAjz*) 



(9) 



where ^^(-2*) = Jq* dz' / H{z') is the comoving angular-diameter distance. One can also define 
the "dilation scale" as introduced in 



Dv{z) := [d\{z)z/H{z)] 



1/3 



(10) 



The 6dF Galaxy Survey and more recently, the WiggleZ team [27j meausred this quantity at 

30 also measured i^^i^^ 



0.44, z = 0.60 and z = 0.73. Percival et al. 



obtain the combined measurement of BAO/CMB for the quantity 



Dv{z) 
Qt 

Dv(zBAo) 



at z 



0.106, z 

z = 0.35. Combining these measurements with the WMAP-7 measurment of la [26 



0.2 and 

one can 



This has been 



given in [28|. From this one can obtain 



Xbao/cmb 



x*c-'x, 



where 



X 



/ -^4^. - 30.95 \ 
^^(0.106) 

DviO.2) ' 
dAjz^) 
Dv{0.35) 

dAjz*) 
DviOAA) 
dAjz*) 

DviO.6) 
dA{z^) 
\ D^(0.73) 



17.55 



10.11 



8.44 



-6.69 



5.45 



(11) 



(12) 
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Figure 1: Left Figure: The la and 2a contour regions in the ilmo — ^so plane. The "dot" represents the best fit 
value. Right Figure: The recosntructed Wde behaviour as a function of redshift. The dahsed line is the best fit 
behaviour. The light grey region is for la confidence level and the dark grey region is for 2a confidence level 



and 



C 



V 



0.48435 
-0.101383 
-0.164945 
-0.0305703 
-0.097874 
-0.106738 



-0.101383 
3.2882 

-2.45497 
-0.0787898 
-0.252254 

-0.2751 



-0.164945 
-2.45497 
9.55916 
-0.128187 
-0.410404 
-0.447574 



-0.0305703 
-0.0787898 
-0.128187 

2.78728 
-2.75632 

1.16437 



-0.097874 
-0.252254 
-0.410404 
-2.75632 
14.9245 
-7.32441 



-0.106738 \ 
-0.2751 
-0.447574 
1.16437 
-7.32441 
14.5022 / 



(13) 

is the inverse covariance matrix derived by using the above results together with the correlation 
coefficients r = 0.337, r = 0.369 and r = 0.438 calculated for the rJDy pair of measurements 
at 2 = (0.2,0.35), z = (0.44,0.6) and z = (0.6,0.73), respectively 



We further use the measurement of CMB anisotropy by WMAP-7 observations [5|]. Using the 
publicly available code CAMB [32,], we calculate the expected CMB anisotrpy spectrum for our 
model. To incorporate the dark energy perturbation, we use publicly available PPF Module [33 
for CAMB. With this we fit our three model parameters ^Imo, ^so and v^. The other cosmological 
parameters are equated to their best fit values as derived from the WMAP-7 data. We compare 
our theoretical results with the data and do the likelihood analysis using the likelihood code 
provided by the WMAP team [sj. 

With these four set of observational data, we constrain our model parameters and also recon- 
struct the dark energy equation of state Wde- The results are shown in Figures (1), (2) and (3) 
where we have set equal to 0, 0.1 and 0.17 respectively which corresponds to Wg = —1/3, —0.27 
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Figure 2: Left Figure: The la and 2u contour regions in the fimo — ^^so plane. The "dot" represents the best fit 
value. Right Figure: The recosntructed Wde behaviour as a function of redshift. The dahsed line is the best fit 
behaviour. The light grey region is for la confidence level and the dark grey region is for 2f7 confidence level 



and —0.22 respectively. With increasing value of , the allowed range for VLso decreases, whereas 
the constraint on VLmo does not change appreciably. Also the contribution from the string net- 
work is always sufficiently small than the rest of the energy density, ensuring that the string 
network never dominates the energy density of the universe. The small contribution from fi^o 
results the Wde to be very close to w = — 1 at present which is evident in the reconstructed be- 
haviour for Wde- As we go past, the equation of state for the dark energy can change appreciably 
from the = — 1 behaviour. In fact for higher values of ff, it is possible that Wde can also scale 
the matter component {wde = 0) at 2a confidence level. 

To summarise, the current observational data is fully consistent with the presence of a small 
but non zero contribution from cosmic string networks. This small contribution is enough to 
produce appreciable change in the behaviour of dark energy equation of state Wde at earlier times 
{z > 0) although one can not distinguish it from w = —1 at the present day. 



Next we calculate the Bayesian Evidence for our model. This is defined as [34 

E = J C{9)P{9)d9, (14) 

where 6 represents the set of model parameters, C represents the Likelihood function and P{0) 
is the prior probability distribution for parameter 9. According to Jeffrey's interpretation jssj, if 
Alnii^ between 1 and 2.5, it is a significant evidence for a model with higher E while the same 
between 2.5 and 5 is a strong to very strong evidence. If AlnE is more than 5, the model with 
higher E is decisively favoured. 
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Figure 3: Left Figure: The la and 2a contour regions in the ilmo — ^so plane. The "dot" represents the best fit 
value. Right Figure: The recosntructed Wde behaviour as a function of redshift. The dahsed line is the best fit 
behaviour. The light grey region is for la confidence level and the dark grey region is for 2a confidence level. 



We calculate the evidence for the simple cosmic string network model for which f ^ = and 
Ws = —1/3. With such choice, we have two model parameters, flmo and Qso- We assume uniform 
prior to these parameters. With this the AlnE for our model and a simple ACDM model is 0.97 
which shows that evidence-wise our model is equally viable as the concordance ACDM model. 

To conclude, we propose a model that can explain the deviation from the concordance ACDM 
behaviour. The majority of the approaches for such a model assume the A term to be zero and 
introduce dynamical fields that can give rise to late time acceleration of the universe. All these 



models are phenomenological in nature and lacks theoretical understanding ( See [36[ for recent 
attempt to build a quintessence model in string theory). They not only inherit the problems 
present in the concordance ACDM model, but also burden us with a number of extra parameters 
which also have to be fine tuned. This motivates us to look for alternative approaches. We 
propose that the presence of frustrated network of cosmic strings together with the already 
existing A can mimick a suitable dark energy model which is consistent with the observational 
data. Although this model shares all the difficulties that are present in a ACDM model, but 
the inclusion of cosmic string network which can cause the observed deviation from the ACDM 
behaviour can be justified as we expect such defect networks to be formed during the phase 
transition in the early universe. 

We put constraint on the model parameters using the observational data and the recon- 
structed equation of state for the dark energy, Wde, shows that, although at present, it is indis- 
tinguishable from w = —1, as one goes to the past, its deviates sufficiently from this behaviour. 
The interesting fact is that to get such large deviation, one needs a very small contribution from 
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the cosmic string network, preventing it to dominate the energy budget of the universe. 

Finally by using Bayesian Evidence, we show that statistically our model and ACDM model 
are equally favored. 
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